Cluster Dynamical Mean-Field Theory of the density-driven Mott transition in the 

one-dimensional Hubbard model 
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The one-dimensional Hubbard model is investigated by means of two different cluster schemes 
suited to introduce short-range spatial correlations beyond the single-site Dynamical Mean-Field 
Theory, namely the Cellular Dynamical Mean-Field Theory, which does not impose the lattice 
symmetries, and its periodized version in which translational symmetry is recovered. It is shown 
that both cluster schemes are able to describe with extreme accuracy the evolution of the density 
as a function of the chemical potential from the Mott insulator to the metallic state. Using exact 
diagonalization to solve the cluster impurity model, we discuss the role of the truncation of the 
Hilbert space of the bath, and propose an algorithm that gives higher weights to the low frequency 
hybridization matrix elements and improves the speed of the convergence of the algorithm. 

PACS numbers: 71.10.-w, 71.27.+a, 75.20.Hr, 75.10.Lp 



Strongly correlated electronic systems and the mod- 
els describing them represent a formidable challenge for 
theorists. An important advance in this field has been 
achieved through the development of the Dynamical 
Mean-Field Theory (DMFT)ji in its single-site version. 
DMFT is a non perturbative approach which fully re- 
tains local quantum dynamics, but simplifies the spatial 
dependence of the correlation functions to make the prob- 
lem tractable. The lattice problem is therefore mapped 
onto a dynamical local problem, and consequently onto 
an impurity model subject to a self-consistency relation. 
The DMFT proved its power providing the first unified 
scenario for the longstanding problem of the Mott transi- 
tion in the Hubbard model, and completely characteriz- 
ing the peculiar criticality associated with this transition. 
Numerous theoretical predictions of this approach have 
been verified experimentally^ The combination of single- 
site DMFT with electronic structure methods has given 
new insights into the physical properties of many corre- 
lated materials. The main limitation of DMFT is the 
neglect of spatial correlation in the one-electron spec- 
tra, which makes impossible to treat phases with definite 
spatial ordering such as d-wave superconductivity. Non 
local effects, such as short-range order, can also be very 
important in describing the metal-insulator transition in 
materials such as Ti2C>3. 3 

There is therefore a strong motivation for developing 
extensions of single-site DMFT. Several schemes have 
been proposed, and they introduce short-range corre- 
lations by replacing the single impurity model with a 
cluster-impurity, containing N c sites with a given spatial 
arrangement. The Dynamical Cluster Approximation 4 
consists in a systematic inclusion of a few lattice mo- 
menta, which would correspond to the lattice momenta 
of an A c -site cluster. From the impurity model point 



of view, this approach requires periodic boundary condi- 
tions on the cluster^ A different approach uses a contin- 
uous interpolation of the sclf-cncrgy in momentum space 
in the self consistency condition** An alternative clus- 
ter method, the Cellular Dynamical Mean -Field Theory 
(CDMFT), has been proposed in Refs. I(l7l closely follow- 
ing the spirit of the DMFT. In this approach, the cluster 
has open boundary conditions, and a 'local' theory for the 
cluster degrees of freedom is obtained through the cavity 
method. The basic approximation is to assume that the 
dynamical field experienced by the cluster is Gaussian. 
A modification of the CDMFT scheme, the Periodized- 
CDMFT (PCDMFT) £ enforces the lattice periodicity at 
every iteration. Both PCMDFT and CDMFT have been 
proved to give causal spectral functions. 8 

In this work we focus on the CDMFT and PCDMFT 
and study their performance in the one-dimensional Hub- 
bard model, which is the worst case scenario for every 
mean-field theory. The ability of CDMFT to capture the 
physics of the Mott insulating state in one dimension has 
been already demonstrated in Ref. |3, where the half- 
filled Hubbard chain was correctly found to be an insula- 
tor for every value of U, with a Mott-Hubbard gap which 
closely follows the Bethe Ansatz (BA) exact solution^ 
Here we extend the analysis of Ref. to the doped metal- 
lic state, and compare CDMFT and PCDMFT to the 
BA results for the whole density range and different val- 
ues of the interaction strength. This is a much more 
stringent test of the method, since metallic states, with 
gapless excitations, are more difficult to describe with a 
reduced number of sites than insulators. One goal of the 
investigation is to compare the performance of CDMFT 
and PCDMFT, with an eye to future applications. As in 
Ref. Oj, we use exact diagonalization to solve the cluster- 
impurity model. A non trivial issue is how to discretize 
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the bath to implement the solution of the cluster impurity 
model by exact diagonalization. This problem was con- 
sidered in Refs llllll2] in the context of single-site DMFT, 
and, more recently, by Potthoff in the context of clus- 
ter methods using a functional technique,— Here we find 
that a simple modification of the notion of distance in the 
algorithm of Ref. ^] improves dramatically the quality 
and the convergence of the results. The paper is orga- 
nized as follows: Sec. II introduces the model and our 
theoretical tools, Sec. Ill presents our results and Sec. 
IV is dedicated to the conclusions. 



I. MODEL AND METHOD 



We consider the one-dimensional Hubbard model 



n^n il - (i rij 



(1) 

where Ci jCT (c icr ) are destruction (creation) operators for 

electrons of spin a, rii a — c\ a ci a is the density of tr-spin 
electrons, t is the hopping amplitude, U is the on-site 
repulsion and /i the chemical potential. 

Within cluster DMFT methods an effective action for 
the cluster degrees of freedom is defined as 



S eff = f drdr'^2cl a (r)g- u Ur-r')c V(T (r') + 

drJ^Un^(T) niil (T), (2) 



l«=l 

where Qz^„ is a Weiss dynamical field and [i, v = 1,...N C 
are indices of sites in the cluster. The action for the clus- 
ter degrees of freedom can be derived, in close analogy 
with DMFT, by means of the cavity construction, i.e., by 
integrating out the degrees of freedom of all the fermions 
except the ones in the cluster, and assuming that the dy- 
namical field due to the other sites is Gaussian^ The clus- 
ter has therefore naturally open boundary conditions, as 
opposed to the periodic boundary conditions of the Dy- 
namical Cluster Approximation^, and breaks the trans- 
lational symmetry. In practical implementations of this 
approach, the initial dynamical Weiss field is specified by 
the choice of a given Q^ a . The effective action is then 
solved with some technique, providing the cluster Green's 
function G Ai;y(T (r) = — (T r c Alcr (r)cJ, (T ). Consequently, the 
cluster self-energy is obtained (dropping the spin index 
to simplify the notations) as 



(3) 



Just like the single-site DMFT, a cluster-DMFT 
scheme is completed by self-consistence relation which 
allows to determine a new Weiss field through the knowl- 
edge of the Green's function. The two methods we com- 
pare in this paper differ in the way this self-consistency 
is implemented. 



A. CDMFT method 

Within CDMFT the cluster self-energy © is used to 
compute the N c x N c matrix of "local" Green's function 
for the cluster, given by 



ir/N c 



1 



dk 



(4) 

where the hat labels N c x N c matrices, the momentum- 
integral extends on the reduced Brillouin zone associated 
to the 7V c -site cluster, tk is the Fourier transform of the 
cluster hopping term^ and 1 is the N c x N c unit matrix. 
Gi oc (iuj n ) and Yi c {iu> n ) are then used to obtain a new 
Weiss field 

(G new )^(i"n) = E^n) + (G lo %l^n), (5) 

which determines the new effective action , from which 
a new cluster Green's function G^ u (ioj n ) is obtained and 
the procedure is iterated until convergence. As noticed 
above, in this approach the cluster is naturally taken with 
open boundary conditions, and breaks the full lattice 
translational invariance. As a consequence, the differ- 
ent sites in the cluster are not equivalent, and the local 
Green's functions on different sites may be different. 



B. PCDMFT method 

The lack of translational invariance of the CDMFT 
may be seen as a limitation of this method. The main 
idea of the PCDMFT approach is therefore to reintro- 
duce the translational invariance in the CDMFT scheme 
without imposing periodic boundary conditions on the 
cluster, as it is done within the Dynamical Cluster 
Approximation. 4 The first step of PCDMFT is to com- 
pute a lattice self-energy with full lattice translational 
invariance by periodizing the cluster self-energy defined 
in Q 



E(k,iu n ) = i 



//// 



E ( u,(«*>n)e 



(6) 



where i? M is the position vector of site fi. Notice that 
the derivation of the lattice self-energy is not unique, 
and alternative estimators have been proposed in Ref. 0- 
The lattice self-energy naturally defines a translationally 
invariant lattice Green's function 



G(k,iuj n ) = 



iuj n + fl-tk — E(fc, iuj n ) ' 



(7) 



where tk is the Fourier transform of the hopping part of 
the Hamiltonian. The key of PCDMFT is to compute the 
local Green's function for the cluster degrees of freedom 
by projecting back the translationally invariant Green's 
function J7J) onto the cluster sites. The PCDMFT ap- 
proximation for the local Green's function of the cluster 
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is therefore given by 



G l °Z(iu n )= I e- ikR »G(k,ioj n )e ikR ^. (S) 



Finally, exactly as in CDMFT, Eq. JSJ) is used to ob- 
tain a new Weiss field and continue the iterative proce- 
dure. We emphasize again that the only difference be- 
tween CDMFT and PCDMFT is in the way the local 
Green's function is computed from the cluster self-energy, 
i.e., in the difference between Eq. @ and Eq. JgJ. The 
latter equation partially reintroduces the full lattice sym- 
metry at each iteration. Finally, it must be noticed that, 
as shown in Ref. H both the CDMFT and PCDMFT 
methods are causal, i.e., they produce self-energies with 
a negative imaginary part. 



C. Exact Diagonalization 

Besides the self-consistency defined by the above equa- 
tions, the CDMFT requires a solution of the cluster ac- 
tion (J2J and the evaluation of the Green's function. For 
practical purposes, it is useful to resort to a Hamiltonian 
formulation, where the quantum fluctuations on the lat- 
tice are realized by hybridization with a conduction bath. 
This leads to a cluster-impurity Hamiltonian of the form 

Hac/ = E » v c t<y C ^ + U X! n ^ n ^ + 

£ k at a a k<J + Y (Vfe 1(U o^.c M<7 + h.c), (9) 

kcr kfia 

where the indices fj,, v = 1, . . . , N c label the cluster sites, 
contains the hopping matrix elements inside the 
cluster and the chemical potential term. The auxiliary 
fcrmionic degrees of freedom a ka describe the bath. In- 
tegrating out the bath degrees of freedom a k a we obtain 
an action of the form with 



iu> n - £ k ' 



(10) 



to create degenerate levels. In particular, we have found 
that each energy level is typically two-fold degenerate 
for the two-site cluster N c = 2. This result suggested 
us to simplify the parameterization assuming that each 
level in the bath has such a degeneracy. This simplified 
parameterization of the bath, which involves fewer pa- 
rameters, typically makes the convergence of (P)CDMFT 
much faster (in terms of the number of iterations needed 
to reach a given threshold) and allows to easily reach ex- 
tremely accurate solution which respect the symmetries 
of the problem. This observation also suggests that some 
values of Nf, should provide worse results than others for 
the only reason that they could not allow for the de- 
sired degeneracy. As an example, odd number of bath 
sites are not the best choice for the two-site cluster. For 
the N c = 3 case, we still find that doubly degenerate 
energy levels spontaneously develop in the iterative pro- 
cess, even if the practical advantage of the degeneracy is 
much smaller than for N c — 2. Further investigation of 
larger N c clusters is needed to fully understand the origin 
of the degeneracy. 

We emphasize that the above discussed truncation of 
the bath is the only approximation introduced in the ex- 
act diagonalization approach to (C)DMFT. At each iter- 
ation, the Weiss field given by Eq. JSJ must be projected 
on the functional space defined by i|10|) . In practice, the 
parameters e k and V ki which minimize a difference A be- 
tween and (|1U|I are used in the next iterations. What 
is the optimal criteria for the choice of these poles has 
been discussed in the literature, where different crite- 
ria have been used. In this work we describe a simple 
modification of the criteria proposed in Ref. ^] which 
substantially improves the convergence of the algorithm. 
We see that a proper definition of the difference is cru- 
cial as far as the convergence with increasing bath size is 
concerned. The substantial independence of the results 
on Nb (beyond a minimum value) observed in Ref. is 
in fact attributed to a definition of functional distance 
which weighted too heavily the high frequency region, 
and we propose an alternative definition of distance that 
eliminates this problem. 



Spv being the Kronecker delta. 

As discussed in Refs. HlHT3 and more recently in Ref. 
fill, to implement the DMFT equations using an exact 
diagonalization solver, it is necessary to discretize the 
problem by parameterizing the hybridization function in 
terms of a finite number of poles. Hence the hybridization 
function (last term in Eq. I|10|)l is approximated by 



(11) 



Here the sums run over a discrete set of fcj = 1, . . . Nb- 
Notice that in general several e ki can be identical, i.e., 
the eigenvalues of the bath can be degenerate. Our nu- 
merical study shows that the self-consistence loop tends 



II. RESULTS 

The one-dimensional Hubbard model represents an 
ideal benchmark for our cluster methods. On the one 
side we can compare with an exact solution, while on 
the other, we know that the single-site DMFT is not 
able to reproduce the Mott insulating state for arbitrarily 
small value of U and the divergence of the compressibility 
K = dn/dfi approaching the Mott transition^ In Ref. 
the ability of CDMFT to reproduce the insulating state 
and the dependence of the Mott-Hubbard gap on U has 
been proven. Here we focus on the behavior of the den- 
sity n as a function of the chemical potential [i and on 
the density-driven Mott transition. 

The n-fi curves for U/t — 1 and U/t = 4 are reported 
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in Figs. Hand|2 respectively. The BA results are com- 
pared with single-site DMFT and (P)CDMFT for N c = 2 
and Nj, = 8. In the weakly correlated case (U/t = 1) 
the difference between the different methods are hardly 
noticeable in the plots, but it must be noted that the 
DMFT gives a metallic solution at half-filling as opposed 
to the insulating exact solution which both CDMFT and 
PCDMFT correctly reproduce already for N c — 2. The 
more correlated U/t — 4 case emphasizes the great im- 
provement on DMFT brought by CDMFT. Besides the 
overall extremely good agreement with BA, CDMFT is in 
fact able to closely follow the compressibility divergence 
close to the Mott transition, which is missed by DMFT 
(upper-left inset of Fig. EJ. PCDMFT is instead unable 
to capture this divergence, and gives a finite compress- 
ibility, but it becomes more and more accurate in the 
extremely doped region (n < 0.5) due to the explicit re- 
inforcement of lattice periodicity (lower- right inset of Fig. 
|2J) • The comparison between the two methods, as well as 
the dependence of the results on the number of bath sites 
is made more quantitative in Fig. |3 where the deviation 
from the exact density is shown for different methods 
and bath sizes for U/t — 4. Not surprisingly, for both 
techniques, the results improve as Nb is increased. In the 
region close to the Mott-Hubbard transition the CDMFT 
results nicely approach the exact solution with k — ► oo 
as Nb is increased, while for the PCDMFT the improve- 
ment with bath size is not able to reproduce this feature. 
Deep in the metallic phase both techniques become ex- 
tremely accurate, but the PCDMFT is typically better 
for n < 0.5. As an example for N c = 2 and Nb = 8, the 
PCDMFT error on the density for fi = -0.5 (n ~ 0.555) 
is 6 x 10~ 4 , while for CDMFT it is 1 x 10~ 3 (See the 
lower-right inset of Fig. 0) . In both cases Nb small as 6 
or 8 is sufficient to achieve remarkable accuracy. We also 
show results for N c = 3 and Nb — 6. In this case, even if 
the cluster size is increased with respect to N c = 2, the 
method has some slight problems due to the odd num- 
ber of sites, which result in a final critical value of U/t 
for the Mott transition (even if extremely reduced with 
respect to DMFT) and a worse agreement with exact re- 
sults close to the Mott transition. However, moving to 
the metallic state, the results rapidly approach the exact 
solution as well as the N c = 2 ones. 

We notice that the expected improvement of results 
with bath size that we have found was not observed in 
Ref. 0- We attribute this difference to the more ef- 
ficient way to discretize the bath we use in this work. 
More precisely, we use a different definition of distance 
between the continuous and the discretized Weiss fields 
that gives more weight to the low Matsubara frequencies 
which arc important to determine low-energy physics, 
namely A = ^ v \Q ris — Q\/uj v - We have checked that, 
using non- weighted differences and a large number of fre- 
quencies, as it was done in Ref. 0, the low-energy part of 
the spectrum has extremely small weight, and it is there- 
fore almost irrelevant. In this case the inclusion of further 
bath sites, which would improve the low energy part of 



the spectrum, is therefore basically irrelevant. We no- 
tice in passing that the use of the weighted distance also 
favored an easier convergence of the iterations. 

These considerations are also supported by the analy- 
sis of dynamical quantities. As in Ref. 0, we consider the 
imaginary part of the on-site Green's function G\\ and 
the real part of the nearest-neighbor Green's function Gn 
on the Matsubara axis, which are plotted for U/t =1 and 
7 and for n = 1. single-site DMFT, CDMFT and the 
PCDMFT for N c = 2 and N b = 8 are compared with 
the results of Density-Matrix Renormalization Group 
(DMRG), a numerical approach which is known to pro- 
vide basically exact results for one-dimensional systems 
(for more details on our calculation of dynamical prop- 
erties with DMRG, see Ref. 0). The agreement of both 
CDMFT and PCDMFT with the virtually exact DMRG 
results is extremely good for U/t = 1, where the pre- 
vious CDMFT results of Ref. were quite inaccurate. 
Interestingly, the single-site DMFT completely fails in 
the description of dynamical properties even if the n-fj, 
curve shown in Fig. ^is close to the exact solution. In the 
strong coupling case U/t = 7, CDMFT closely follows the 
DMRG, while PCDMFT is a poorer approximation, but 
still substantially better than the single-site DMFT. We 
also experienced some problem in obtaining converged 
PCDMFT solutions for U/t > 7. 



III. CONCLUSIONS 

We have tested the ability of CDMFT and PCDMFT 
to capture the physics of the one-dimensional Hubbard 
model both in the metallic and the insulating phase. 
While we expected cluster methods to give reasonable 
results for short-distance high-energy physical quantities, 
we found that both methods are remarkably accurate for 
basically all the metallic densities, where the low-energy 
physics is important. The CDMFT method naturally 
generalizes DMFT to include short-range dynamical cor- 
relations within an open cluster of N c sites, which breaks 
the full lattice translational symmetry. PCDMFT par- 
tially restores the translational symmetry without im- 
posing periodic boundary conditions to the cluster. The 
CDMFT, already for N c — 2, is able to reproduce the 
divergent compressibility of the one-dimensional system 
when the Mott transition is approached, providing a 
qualitative change with respect to single-site DMFT. The 
PCDMFT, due to the restored translational invariance, 
gives slightly better results than CDMFT deep in the 
metallic regime (far from the Mott transition), where 
correlation effects are smaller, but fails in reproducing 
correctly the physics close to the Mott transition. The re- 
markably good results of CDMFT (either from the metal- 
lic to the insulating state) and of PCDMFT (in the metal- 
lic state only) are also due to a technical improvement in 
the truncation procedure of the bath hybridization func- 
tion of the impurity problem. A systematic improvement 
of the accuracy of the results by increasing the bath size 
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is observed. We also find that some degeneracy of the 
bath energy levels automatically develops if not explic- 
itly imposed. Explicit enforcement of degenerate levels 
with the correct degeneracy strongly favors convergence 
and accuracy of the solutions. These methodological ad- 
vances will be useful in further studies of more complex 
systems. 
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FIG. 1: Density n as a function of /i for U/t = 1 N c = 2,Nb = 
8 within single-site DMFT, CDMFT and PCDMFT compared 
with the exact solution. 




FIG. 2: Density n as a function of fj, for U/t = 4, N c — 2,N b = 
8 . The upper-left inset show a close up of the region close to 
half-filling, and the lower-right inset shows a region deep in 
the metallic regime. 
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FIG. 3: Deviations between the CDMFT and PCDMFT den- 
sities n with respect to the exact solution (BA) for U/t — 4. 
In the upper/lower panel we show CDMFT/PCDMFT for 
different bath and cluster sizes (see legend). 




FIG. 4: Imaginary part of the local Green's function Gu and 
real part of nearest-neighbors Green's functions G12 for n = 1. 
Single-site DMFT and two-site CDMFT and PCDMFT are 
compared with DMRG for U/t = 1 (left) and U/t = 7 (right). 



